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Exergoeconomics  is  also  called  thermoeconomics,  and  thermoeconomic  analysis  methodologies  combine 
economic  and  thermodynamic  analysis  by  applying  the  cost  concept  to  exergy  which  accounts  for  the 
quality  of  energy.  The  main  concept  of  thermoeconomics  is  the  exergetic  cost  and  it  deals  with  cost 
accounting  methods.  This  paper  is  a  review  on  the  exergoeconomic  analysis  and  optimization  of  combined 
heat  and  power  production  (CHPP).  A  brief  historical  overview  on  the  exergoeconomics  analysis  and 
optimization  is  given.  The  concept  of  exergetic  cost  and  cost  accounting  methods  are  discussed.  An 
application  of  relevant  formulation  is  given  using  a  diesel  engine  powered  cogeneration  system  as  an 
example.  Main  thermoeconomic  methodologies  available  in  literature  are  described  and  their  advantages 
and  disadvantages  with  respect  to  one  another  are  compared  and  discussed  through  a  well-known 
problem,  namely  CGAM.  Important  studies  on  thermoeconomic  analysis  and  optimization  of  combined 
heat  and  power  production  are  listed  based  on  the  methodology  used  and  the  type  of  system  considered. 
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Nomenclature 

c 

cost  per  unit  of  exergy  ($/GJ) 

C 

cost  rate  associated  with  exergy  ($/h) 

E 

exergy  rate  (kW) 

^dest 

rate  of  exergy  destruction  (kW) 

^heat 

rate  of  exergy  transfer  by  heat  (kW) 

/ 

exergoeconomic  factor  (%) 

h 

enthalpy  (kj/kg) 

I 

investment  cost  ($) 

L 

Lagrangian 

L 

rate  of  the  Lagrangian 

n 

exergoeconomic  cost  vector  (kW) 

r 

relative  cost  difference  (%) 

S 

entropy  rate  (kW/K) 

Q 

rate  of  heat  transfer  (kW) 

w 

power  (kW) 

X 

vector  of  decision  variables 

y 

vector  of  all  the  inputs  and  outputs 

Y 

constraint  function 

z 

cost  rate  associated  with  the  sum  of  capital 
investment  and  O&M  ($/h) 

Abbreviations 

TADEUS 

thermoeconomic  analysis  for  diagnosis  of  effi¬ 
ciency  reduction  in  systems 

DEPC 

diesel  engine  powered  cogeneration 

TEC 

theory  of  the  exergetic  cost 

TECD 

theory  of  the  exergetic  cost— disaggregation  meth¬ 
odology 

TFA 

thermoeconomical  functional  analysis 

EEA 

exergoeconomical  analysis 

SPECO 

specific  exergetic  cost  analysis  method 

AVCO 

average  exergetic  cost  analysis  method 

MOPSA 

modified  productive  structural  analysis  method 

LIFO 

last  in  first  out  principle 

EFA 

engineering  functional  analysis 

STT 

structural  theory  of  thermoeconomics 

C 

compressor 

AP 

air  preheater 

T 

turbine 

CC 

combustion  chamber 

Ev 

evaporator 

Ec 

economizer 

HRSG 

heat  recovery  steam  generator 

Greek  letters 

l 1 

energy  efficiency 

8 

exergy  efficiency 

P 

capital  recovery  factor  (%) 

X 

vector  of  all  the  Lagrange  multipliers 

X 

number  of  units  and  junctions 

CO 

number  of  units,  junctions  and  branching  points 

Superscripts 

BQ 

steam 

CHE 

chemical 

P 

mechanical 

T 

thermal 

W  work  or  electricity 

Subscripts 
av  average 

D  exergy  destruction 

F  fuel 

i  any  stream 

in  streams  entering  a  component  or  system 

k  any  component 

mar  marginal 

out  streams  leaving  a  component  or  system 

P  product 

tot  total  plant 


1.  Introduction 

Exergoeconomics  (thermoeconomics)  combines  exergy  analy¬ 
sis  with  conventional  cost  analysis  in  order  to  assess  and  improve 
the  performance  of  energy  systems  [1].  The  primary  contribution 
of  an  exergy  analysis  to  the  evaluation  of  an  energy  system  comes 
through  a  thermoeconomic  evaluation  which  considers  not  only 
the  inefficiencies  but  also  the  costs  associated  with  these 
inefficiencies  and  the  investment  expenditures  required  to  reduce 
them.  The  objective  of  a  thermodynamic  optimization  is  to 
minimize  the  thermodynamic  inefficiencies  within  the  system, 
whereas  the  objective  of  a  thermoeconomic  optimization  of  a 
system  is  to  estimate  the  cost-optimal  structure  and  the  cost- 
optimal  values  of  the  thermodynamic  inefficiencies  in  each 
component  [2-4]. 

In  this  paper,  starting  from  earlier  pioneer  works  in  the  area  of 
thermoeconomics,  major  contributions  are  summarized.  By 
detailed  search  with  appropriate  keywords,  a  researcher  can 
easily  reach  to  the  related  papers  in  the  open  literature.  First,  a 
short  historical  overview  on  the  exergoeconomic  analysis  and 
optimization  of  thermal  systems  is  given.  In  the  following  sections, 
exergy  costing  and  cost  accounting  concepts  are  described.  Various 
thermoeconomic  methodologies  are  discussed  and  exergoeco¬ 
nomic  analysis  and  optimization  of  combined  heat  and  power 
production  are  presented. 

2.  Historical  overview 

The  idea  of  coupling  exergy  and  cost  streams  was  first 
discussed  by  Keenan  in  1932  [5].  He  pointed  out  that  the  value 
of  the  steam  and  the  electricity  rests  in  the  “availability”  not  in 
their  energy  [6-8].  In  the  late  1950s,  the  studies  of  second  law 
costing  started  in  two  different  places  independently.  Tribus  and 
Evans  [6]  studied  desalination  processes  by  exergy  analysis, 
which  led  them  to  the  idea  of  exergy  costing  and  its  applications 
to  engineering  economics,  for  which  they  coined  the  word 
“ thermoeconomics ”  [9].  The  concept  of  their  procedure  was  to 
trace  the  flow  of  money,  fuel  cost  and  operation  and  amortized 
capital  cost  through  a  plant,  associating  the  utility  of  each 
stream  with  its  exergy.  In  the  late  1960s,  El-Sayed  connected 
with  Evans  and  Tribus  in  their  research  on  desalination  and  they 
published  one  of  the  important  papers  in  the  subject  in  1970,  in 
which  the  mathematical  foundation  for  thermal  system  opti¬ 
mization  was  given  [10].  Another  study  on  their  idea  was 
performed  by  Reistad  who  applied  the  method  of  El-Sayed  and 
Evans  to  a  simple  power  plant,  comparing  that  approach  with 
conventional  optimization  procedures  [11].  Also  in  1960s, 
Gaggioli  studied  the  optimal  design  of  power  plant  steam 
piping  and  its  insulation  in  his  Ph.D.  thesis.  He  proposed  costing 
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steam  exergy  at  a  value  to  that  of  power  produced,  penalizing 
exergy  destructions  and  losses  for  the  electricity  which  there¬ 
fore  would  not  be  produced  [12]. 

In  Europe,  many  important  works  on  the  second  law  analysis 
methodologies  and  on  “exergy”  itself  were  performed  in  late  1950s 
and  during  1960s  and  1970s  [13-20].  Bergmann  and  Schmidt 
assigned  costs  to  the  exergy  destruction  in  each  component  of  a 
steam  power  plant  and  optimized  feed  water  heaters  [17].  Szargut 
used  exergy  costing  procedure  in  the  analysis  of  a  simple 
cogeneration  plant  by  derivation  of  an  approximate  generalized 
formula  and  introduced  ecological  cost  coefficient  into  the 
literature.  The  use  of  this  coefficient  made  the  determination  of 
cumulative  consumption  of  non-renewable  natural  resources  in  a 
production  process  possible  [18-20].  Some  of  these  works  were  on 
thermoeconomic  analysis  and  collected  in  a  book  in  1985  by  Kotas 
[21  ],  which  is  still  accepted  as  one  of  the  basic  references  in  exergy 
analysis  and  thermoeconomics  of  thermal  systems. 

Although  there  had  been  numerous  studies  and  theoretical 
approaches  in  the  field  of  thermoeconomics,  the  methodological 
and  functional  applications  of  it  to  the  analysis,  design  and 
optimization  of  thermal  systems  comprehensively  did  not  start 
until  1980s.  During  this  decade,  the  interest  and  studies  related  to 
exergoeconomic  analysis  techniques  and  applications  have 
increased  considerably.  Frangopoulos  [22-27]  and  von  Spakovsky 
[26-30]  applied  and  formalized  the  first  autonomous  method 
developed  by  Evans  and  El-Sayed  [10].  Tsatsaronis  introduced 
some  key  concepts  of  thermoeconomics  such  as  Fuel  and  Product 
[31].  Valero  et  al.  published  another  key  paper  on  thermoeco¬ 
nomics  and  presented  the  basic  methodology  related  to  exergy 
based  cost  analysis  and  applications  [32]  and  with  Lozano,  he 
presented  basics  and  several  applications  of  the  theory  of  exergetic 
cost,  a  major  cornerstone  approach  to  the  field  of  thermoeco¬ 
nomics  [33]. 

The  major  contributions  were  done  in  1990s,  to  achieve  a 
greater  standardization  and  formalism  in  the  area  of  thermo¬ 
economic  studies.  The  common  idea  was  to  propose  a  standard  and 
common  mathematical  formulation  for  all  thermoeconomic 
methodologies  employing  thermoeconomic  models  that  can  be 
expressed  by  linear  equations  [34].  At  that  time,  Tsatsaronis 
proposed  the  term  “ exergoeconomics ”  which  was  defined  as  a  part 
of  thermoeconomics.  Since  the  latter  have  been  used  in  a  general 
sense  expressing  the  interaction  between  any  thermodynamic 
variables  and  economics,  he  suggested  that  exergy  based  cost 
accounting  methodologies  should  be  indicated  by  exergoeco¬ 
nomics  [35].  One  of  the  most  interesting  works  was  the  CGAM 
problem  [36-40]  developed  originally  by  Frangopoulos,  Tsatsar¬ 
onis,  Valero  and  von  Spakovsky.  The  CGAM  problem  was  named 
after  the  first  initials  of  the  participating  investigators.  The 
objective  of  the  problem  was  to  show  how  the  methodologies 
were  applied,  what  concepts  were  used  and  which  results  were 
obtained  in  a  simple  and  specific  problem  [36].  The  sequential 
papers  published  in  the  same  issue  had  the  common  definitions  of 
the  physical,  thermodynamic  and  cost  models  plus  the  objective 
function  but  different  procedures  of  thermoeconomic  optimiza¬ 
tion.  The  authors  emphasized  that  the  aim  of  the  presented  work 
was  the  unification  of  thermoeconomic  methodologies,  not  a 
competition  among  them.  Another  important  project  called 
TADEUS  (in  honor  of  Tadeus  Kotas  [21])  was  initiated  in  2001 
[41-43].  The  major  goals  of  that  project  were  to  apply  procedures 
from  different  research  groups  in  thermoeconomic  analysis  to  the 
diagnosis  of  the  energy  system  malfunction  and  inefficiencies  and 
to  establish  the  common  concepts  and  nomenclature  and  compare 
the  results  and  highlight  the  main  characteristics  of  each  approach. 

Since  1980s,  there  have  been  numerous  published  papers  all 
around  the  world  on  exergoeconomic  cost  analysis,  application 
and  optimization  of  thermal  systems.  Good  fractions  of  them  have 


been  published  since  mid-90s  due  to  the  improved  structural 
formalism  of  the  exergoeconomic  methodologies  [44-81].  Exer¬ 
goeconomic  studies  have  been  successfully  applied  to  power 
plants,  combined  heat  and  power  production  or  cogeneration 
facilities  [45-47,50,51,54,57,59-81].  Progressed  innovative  meth¬ 
odologies  of  fuzzy  logic  and  genetic  algorithm  methods  have  also 
been  applied  to  existing  power  plants  and  cogeneration  facilities 
since  late  1990s  [47,59,69]. 


3.  The  cost  of  exergy  and  exergoeconomics 


The  exergy  cost  of  a  mass  or  energy  stream  is  the  amount  of 
exergy  required  to  produce  it.  In  the  case  of  a  cogeneration  plant, 
the  exergy  cost  of  the  net  power  is  the  exergy  provided  by  fuel  to 
generate  the  electrical  power  delivered  to  the  network  by  the 
cogeneration  plant  [33-39].  When  analyzing  the  cost  structure  of  a 
system,  one  can  distinguish  between  average  costs  and  marginal 
costs  [1-3,33].  Average  costs  are  ratios  and  express  the  average 
amount  of  resources  per  unit  of  product,  and  can  be  defined 
mathematically  as  [2,8,21,33] 


Cavd  — 


c 


av,t 


Ei 


where  cav  i  is  the  unit  specific  average  cost  of  the  ith  stream,  Cavd 
and  Ei  represent  the  exergetic  cost  and  exergy  of  the  ith  stream, 
respectively.  Marginal  costs  are  derivation  and  indicate  the 
additional  resources  required  to  generate  one  or  more  units  of 
product  under  specified  conditions.  Mathematically,  it  is  defined  as 
[33] 


Cmar,i 


Since  the  average  cost  is  not  predictive,  it  can  only  be  known 
after  production  process,  when  we  know  how  much  resources 
were  used.  On  the  other  hand,  marginal  costs  are  predictive  in 
nature;  they  can  be  used  to  calculate  additional  fuel  consumption 
when  the  production  is  modified.  Thermoeconomic  optimization 
methods  are  generally  based  on  the  marginal  costs  while  average 
costs  are  used  in  the  exergoeconomic  analysis  of  the  systems  [32- 
34].  In  exergy  costing,  a  cost  is  associated  with  each  exergy  stream. 
Exergy  transfers  by  the  entering  and  exiting  streams  of  matter  and 
by  power  and  heat  transfer  rates  may  be  written,  respectively  as 
[2] 

Q  =  CiEi  (3) 

Ce  —  ceEe  (4) 


CW  -  CyyVV 


Cq  —  cqEq  (6) 

Accordingly,  for  a  kth  component  receiving  heat  transfer  and 
generating  power,  we  may  write  [2] 

^  ^  (CeEe)k  +  Cw  fcVVk  —  CqfcEqk  +  ^  ^  (7) 

e  i 

In  the  case  of  a  single  product  of  the  plant,  the  unit  cost  of  the 
product  can  be  easily  determined  using  Eq.  (7).  In  the  case  of  a 
multi-product  plant  such  as  cogeneration  plants,  the  cost-balance 
equation  is  not  sufficient.  Additional  criteria  are  required  to 
determine  the  relationship  between  the  unit  costs  of  the  different 
products.  This  is  where  exergy  can  be  used  as  a  basis  for  cost 
allocation  of  the  products. 

A  diesel  engine  powered  cogeneration  (DEPC)  plant  (see 
Fig.  1)  is  considered  as  an  illustrative  example  of  the  exergy 
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Fig.  1.  The  Schematic  of  an  actual  diesel  engine  powered  cogeneration  plant.  C:  compressor,  T:  turbine,  WHB:  waste  heat  boiler,  DeSOx:  desulphurization,  AWR:  air-water 
radiator,  IC:  intercooler,  LOC:  lubrication  oil  cooler,  LOT:  lubrication  oil  tank,  SD:  steam  drum,  FWT:  feed  water  tank,  FDT:  fuel  oil  day  tank,  CON:  condenser,  FFM:  fuel 
forwarding  module,  and  P:  pump. 


costing  of  co-generated  products.  When  the  engine  starts,  air  is 
charged  to  the  compressor  of  the  turbocharger  unit.  The 
turbocharger  consists  of  a  turbine  and  a  compressor,  which 
are  mounted  on  a  common  shaft.  The  shaft  power  obtained  by 
the  operation  of  diesel  engine  is  transferred  to  the  generator  for 
electricity  production.  The  exhaust  gases  leaving  the  engine  flow 
through  the  turbine  of  the  turbocharger  unit  to  produce  the 
necessary  shaft  work  for  the  compressor.  The  air  leaving  the 
compressor  is  cooled  by  water  in  an  intercooler  before  air  enters 
the  engine  cylinders.  The  exhaust  gases  leaving  the  turbine  enter 
the  exhaust  gas  boiler  unit  to  transfer  heat  to  the  feed  water  to 
produce  steam  for  manufacturing  facilities  in  the  factory  and  for 
preheating  of  streams  in  the  auxiliary  equipments  such  as  fuel 
forwarding  module  (FFM)  and  fuel  oil  in  daily  usage  tank  (FDT). 
The  exhaust  gas  leaving  the  boiler  is  sent  to  the  DeSOx 
(desulphurization)  unit  in  which  the  SOx  emission  is  lowered 
to  the  acceptable  legal  values.  After  the  DeSOx  unit,  the  exhaust 
gases  are  released  to  the  atmosphere.  In  Table  1,  exergetic  cost 
rate  balances  and  related  auxiliary  equations  for  each  subsystem 
of  the  DEPC  plant  are  presented  [80,81]. 

The  explicit  determinations  of  auxiliary  equations  based  on 
exergy  criteria  were  presented  in  previous  works  [7,8,12]. 
These  methods  are  known  as  the  equality,  the  extraction  and 
the  by-product  methods.  All  these  methods  require  a 
judgment  regarding  the  purpose  of  each  unit  of  the  plant  and 
distribute  the  costs  proportionally  to  exergy  contents  of  the 
flow.  These  three  methods  are  applied  to  the  turbine  of  the 
turbocharger  unit  of  the  DEPC  plant  (see  Fig.  1)  as  an  illustrative 
example.  The  exergetic  cost-balance  equation  of  the  turbine 
from  Table  1  is 

Cq  —  C7  +Zturbine  =  Cwturbine  (8a) 


an  assumption  must  be  made  regarding  the  allocation  of  the  cost  to 
the  two  outputs:  shaft  power  Ew  and  the  exit  stream  E7. 


(i)  In  the  equality  method,  the  generation  of  the  two  products  is 
considered  to  have  the  same  priority,  so  the  cost  of  high 
temperature  exhaust  gas  of  the  turbine  is  charged  to  the 
production  of  both  outputs  (products)  proportionally  to  their 
exergies  as 


c7  =  Cw 


CqEq  +  ZjurbinE 
£7  +  Ew 


(ii)  In  the  extraction  method,  it  is  considered  that  the  purpose  of  the 
turbine  is  to  generate  shaft  power  and  thus  the  entire  cost  of  the 
turbine  must  be  charged  against  it.  This  results  in  assuming  the 
unit  exergetic  cost  of  the  exhaust  gas  entering  the  turbine  to  be 
the  same  as  that  of  the  exhaust  gas  leaving  it.  Thus, 

c7  =  c6  (10) 


Cq  (Eg  —  E 7)  +  ZjURBINE 

CW  =  - ■ -  [  1  1  J 

Ew 

(iii)  In  the  by-product  method,  the  cost  of  one  of  the  outputs  is 
assumed  to  be  known.  It  is  assumed  that  the  generation  of 
process  stream  is  essential  even  if  no  shaft  power  would  be 
generated.  Therefore,  the  exergy  of  the  exhaust  gas  is  costed  as 
if  it  were  produced  alone  in  a  diesel  engine  at  the  required 
pressure  and  temperature.  Therefore,  the  cost  of  shaft  power  is 
determined  by  the  cost  balance. 

Table  1  shows  the  exergetic  cost  rate  balances  and  correspond¬ 
ing  auxiliary  equations  for  each  subsystem  of  the  DEPC  plant  given 
in  Fig.  1  [80,81]  using  the  extraction  method. 


or 

c6E6  +  ZTURBine  =  c7^7  +  cw^w  (8b) 

Assuming  the  exergetic  cost  of  the  inlet  stream  c6  is  known,  the 
cost-balance  equation  has  two  unknowns,  namely,  c7  and  cw.  Thus, 


3.2.  Exergetic  cost  formation  process  and  cost  accounting 
methodologies 

3  A  A.  Fuel-product  model 

Exergetic  cost  accounting  provides  a  wide  and  clear  vision  of  the 
usage  and  degradation  of  energy  and  in  consequence  of  natural 
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Table  1 

Exergetic  cost  rate  balances  and  corresponding  auxiliary  equations  for  each  subsystem  of  the  DEPC  plant.  State  numbers  refer  to  Fig.  1 . 


Component 

Exergetic  cost  rate  balance  equations  for  the  components  of  DEPC 

Auxiliary  equations 

Compressor 

^wCOMP  +  ZCOMP  =  C2  —  Ci 

Cl  =  0  (El 

=  0) 

Intercooler 

(C20  —  C21)  +  (C23  -  C24)  +  ZIC  =  C3  -  C2 

C21  -  C20  

C24  -  C23 

E21  -  E20 

E24  -  E23 

Lubrication  oil  cooler 

(C24  -  C25)  +  Zioc  =  C28  -  C27 

C28  _  C27 
E28  E27 

c  24  =  C25 

C22  _  C21 

Diesel  engine 

Q  +  C5  +  (C21  -  C22)  +  (C26  -  C27)  +  ZDE  =  Cq  +  CW-ELEctric 

E22  E21 

C26  C27 

E26  E27 

Turbine 

Cq  —  C7  +  ZturbINE  =  CwTURBINE 

C6  =  C7 

Waste  heat  boiler 

(C7  —  Cg)  +  ZWhb  =  Cl 3  -  C12 

C12  _  Ci  3 

Ei  2  E13 

C7  =  ^8 

DeSOx  unit 

(Cg  —  Cg)  +  ZDeSO*  =  0 

Cg  =  0 

Fuel  oil  day  tank 

(Cl 5  -  Cl 6)  +ZEDT  =  C4  -  Co.FO 

Ci 5  _  Ci6 

El  5  Ei  6 

Fuel  forwarding  module 

(Ci 7  -  Cig)  +ZEEM  =  C5  -  Q 

Cl  7  _  Cl  8  _ 
El 7  ElS 

Cl  6 

El  6 

Condenser 

(C29  -  C30)  +Zcon  =  C19  -  Cl 8 

Cl 8  _  Cig 

El  8  Eig 

C29  =  0 

Pump 

^Wpyjyjp  ~ 1 —  ^ P  ^11  C,0 

Cio  _  C11 
E10  En 

DEPC 

Cl  +  C5  +  C10  +  ZDEpc  =  Cl  4  +  Cl  5  +  Cg  +  CW-ELEctric 

Ci  4  _  Ci  5  _ 
El  4  E15 

C12 

Ei  2 

resources  [34].  Cost  accounting  methodologies  are  based  on  cost 
allocation  rules.  In  order  to  identify  the  process  of  cost  formation, 
general  relationships  are  required  that  relate  the  overall  efficiency  of 
the  plant  and  the  cost  of  the  products  with  the  efficiency  and  the 
irreversibilities  of  each  individual  component  that  form  the  system 
[33,34,48,61]. 

In  order  to  build  up  the  cost  formation  in  a  system,  productive 
structure  information  of  the  system  must  be  obtained.  This 
structure  consists  of  a  qualitative  and  quantitative  analysis  of  the 
relations  that  link  the  flows  of  exergy  inputs  and  outputs  between 
components  of  an  economic  unit.  Since  all  thermal  plants  have  a 
defined  aim  to  produce  one  or  more  products,  the  quantity  of 
resources  must  also  be  identified  through  mass  or  energy  flows 
which  are  known  as  fuel.  Thus,  each  of  the  components  of  the  plant 
can  be  characterized  by  its  fuel  and  its  product(s).  In  Table  2,  the 


Table  2 

Definitions  of  the  exergies  of  the  fuels  EF  and  the  exergies  of  products  EP  for  the 
components  of  the  DEPC  plant.  State  numbers  refer  to  Fig.  1  [82]. 


Component 

Ef 

EP 

Compressor 

Ec 

E2-E1 

Intercooler 

E2-E3 

(E21  -  E20)  +  (E24  -  E23) 

Lubrication  oil  cooler 

E27  -  E28 

E25  -  E24 

Diesel  engine 

ECf  +  £3  +  £5 

+  E21  +  E26 

E6  +  E22  +  E27  +  E31 

Turbine 

E&  -  E7 

Et 

Waste  heat  boiler 

E7  -  Eg 

E13  Ei 2 

Fuel  oil  day  tank 

E15  -  Ei6 

E4  -  Eo,fo 

Fuel  forwarding  module 

Ei  7  -  Ei  g 

e5  -  E4 

Condenser 

El  8  -  Ei  9 

E30  -  E29 

Pump 

EP 

E11  -  E10 

fuel  and  product  exergy  definitions  of  each  component  in  the  DEPC 
plant  is  given  [82]. 

In  accordance  with  the  fuel-product  model,  the  product  of  one 
component  is  used  as  fuel  for  another  component  or  as  part  of  the 
total  production  of  the  plant.  Thus, 


n 

Ep,i  =  £f,i,0  +  i  =  0,  1 , . . . ,  n 

j= 1 


(12) 


where  EFij  is  the  production  portion  of  the  ith  component  and  is 
the  fuel  of  the  jth  component.  In  Eq.  (12),  the  component  “0”  is 
considered  as  the  environment,  and  then  EFi0  represents  the 
production  portion  of  the  ith  component  which  leads  to  the  final 
product,  coming  from  the  environment  to  the  ith  component.  The 
resources  entering  each  component  can  also  be  expressed  as 


n 

Efj  =  EF, o,i  +  TAn  i  =  0, 1 , . . . ,  n 
j= i 


(13) 


where  EF  O  i  represents  the  external  resources  entering  to  the  plant, 
which  enters  through  ith  component.  The  total  fuel  and  product  of 
the  system  can  be  obtained  as  [34-36,39] 


n 

&F, total  =  ^  ^Ep  Q j  (14) 

j= 1 


n 

Ep,  total  =  VA.W 

j= 1 


(15) 


3.2.2.  The  cost  equations  of  fuel-product  model 

The  proposed  fuel-product  model  can  be  applied  to  exergetic 
cost  structure  formation  of  a  system  and  exergetic  costs  of  product 
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and  fuel  can  be  expressed,  respectively  as 

n 

Cp,i  =  Q,i, o  +  Z^F,ij  i  =  0,1, . . .  ,n 
j= i 


is  similar  to  Frangopoulos’  Thermo  economical  Functional  Analysis 
{ TFA )  [22-25].  The  exergetic  cost  determined  is  distributed  among 
all  other  components  according  to  the  variation  of  entropy  within 
each  one. 


Q v  —  Q,o,/  +  i  —  0, 1 ,  •  •  • ,  n 

i=i 


(17) 


The  proposed  cost  model  of  fuel-product  approximation  is  the 
base  of  the  “Exergetic  Cost  Theory”  [33].  It  is  explained  in  the  next 
section  with  the  cost  assessment  rules  of  the  theory.  The  exergetic 
cost-balance  equations  and  auxiliary  relations  of  the  DEPC  plant 
given  in  Table  1  are  based  on  the  fuel-product  model  and  obtained 
by  using  SPECO  approach  that  is  one  of  the  important  thermo- 
economic  methodologies  described  in  the  next  section. 


4.  Thermoeconomic  methodologies 

Following  El-Sayed  and  Gaggioli  [7,8],  thermoeconomic  meth¬ 
ods  can  be  classified  in  two  groups:  algebraic  methods  and  calculus 
methods.  These  two  basic  categories  are  examined  with  their 
subcategories  and  related  works. 


4.2.  Algebraic  methods 

The  algebraic  methods  use  algebraic  cost-balance  equations 
derived  from  conventional  economic  analysis  and  auxiliary  cost 
equations  for  each  subcomponent  of  any  system  presented  [2]. 
They  are  related  with  the  cost  formation  process  of  the  system  in 
order  to  investigate  the  average  costs. 

4.1.1.  The  theory  of  the  exergetic  cost  (TEC) 

Theory  was  developed  by  Lozano  and  Valero  [33]  and  the 
methodology  presented  in  theory  is  based  on  a  set  of  propositions 
[2].  This  theory  permits  the  introduction  of  a  new  thermodynamic 
concept  called  exergy  cost.  For  given  a  system  whose  limits, 
disaggregation  level,  and  production  aim  of  the  subsystems  have 
been  defined,  exergy  cost  E *  of  a  physical  flow  is  defined  as  the 
amount  of  exergy  needed  to  produce  this  flow  [33].  The  first  step  in 
the  application  of  this  methodology  is  the  division  of  the  system  into 
units,  which  may  adapt  to  a  component  or  a  set  of  components.  A 
single  product  and  fuel  for  each  component  in  the  system  must  be 
defined.  Thus,  a  system  of  equations  can  be  built  with  a  cost-balance 
equation  for  each  unit  (proposition  1 ),  cost  allocation  equations  for 
external  flows  into  the  system  for  which  costs  are  externally  defined 
(proposition  2),  and  losses  for  which  the  cost  is  set  equal  to  zero 
(proposition  3).  However,  this  procedure  is  rarely  enough.  The 
solution  may  be  possible  when  two  following  propositions  are 
considered:  ( 1 )  if  definition  of fuel  of  a  component  includes  a  stream 
that  goes  through  another  component  and  is  used  in  it,  then  the  unit 
cost  of  the  stream  flowing  into  and  out  of  the  component  is  the  same 
(proposition  4a).  (2)  If  the  product  of  a  component  is  composed  of 
two  or  more  streams  then  the  unit  cost  of  those  streams  are  equal 
(proposition  4b).  Using  the  F  (fuel)-P  (product)-L  (loss)  definitions, 
corresponding  matrices  can  be  developed.  Starting  from  these 
matrices  and  using  the  data  from  design  and  operation  of  the  any 
thermal  power  plant,  it  is  possible  to  carry  out  the  exergetic  and 
energetic  analysis  of  the  plant.  The  unit  exergetic  costs  obtained 
from  the  procedure  can  be  used  for  the  optimization  of  any  thermal 
plant.  This  methodology  is  straightforward  in  its  application, 
allowing  for  improvements  and  modifications. 

4.1.2.  The  theory  of  exergetic  cost— disaggregating  methodology 
(TECD  ) 

This  theory  was  proposed  by  Valero  and  Lozano  [33,34] 
searching  the  thermoeconomic  unification  methodology  which 


4.1.3.  Exergoeconomic  analysis  (EEA)  methods 

The  exergoeconomics  approach  was  proposed  by  Tsatsaronis 
and  co-workers  [1-4,35-37,52].  There  are  two  possible  variants  in 
these  methodologies:  specific  cost  and  average  cost.  The  average 
cost  concept  is  very  similar  to  the  application  of  Theory  of  the 
Exergetic  Cost.  In  the  specific-cost  method,  the  cost  of  the  addition  of 
exergy  to  a  material/stream  current  is  determined  and  charged  to 
the  unit  that  makes  use  of  that  exergy.  This  means  that  a 
component  will  obtain  the  exergy  from  a  stream  at  different  costs, 
depending  on  the  components  that  supplied  the  exergy  to  that 
stream.  Due  to  the  purpose  of  the  energy  conversion  plant  (i.e. 
power  plant  or  combined  heat  and  production  plant)  the  division 
between  subcomponents  of  the  main  plant  may  be  made.  If  power 
produced  is  the  primary  product,  cost  of  all  external  irreversi¬ 
bilities  is  loaded  on  the  electricity.  If  heat  is  produced  as  primary 
one,  the  similar  procedure  is  followed  by  the  steam  produced  in  the 
plant.  The  optimization  procedure  of  the  exergoeconomics  is  based 
on  an  iterative  design  improvement  procedure  that  does  not  aim  at 
calculating  the  global  optimum  of  a  predefined  objective  function 
[35],  as  the  conventional  optimization  methods  perform,  but  tries 
to  find  a  good  solution  for  the  overall  system  design.  The  methods 
involve  some  characteristic  parameters  of  exergoeconomic  ana¬ 
lysis  such  as  relative  cost  difference,  exergoeconomic  factor,  and 
exergetic  efficiency.  These  methods  can  be  grouped  in  three 
classes:  (i)  Last-in-First-out  (LIFO)  principle,  (ii)  Specific  Exergy 
Costing/Average  Cost  (SPECO/AVCO)  approach,  and  (iii)  Modified 
Productive  Structure  Analysis  (MOPSA)  approach. 

Last-in-First-out  {LIFO)  principle.  This  approach  was  developed 
by  Tsatsaronis  et  al.  [  1  ].  LIFO  principle  of  cost  accounting  is  used  for 
the  spent  exergy  units  to  calculate  the  cost  of  exergy  supply  to  the 
carrier.  The  method  eliminates  the  need  for  auxiliary  assumptions 
in  the  exergoeconomic  analysis  of  energy  systems  and  improves 
the  fairness  of  the  costing  process.  This  approach  is  based  on  a 
simple  idea,  which  is  directly  linked  between  asset  management 
and  valuation  method:  The  exergy  units  that  are  supplied  last  to 
the  material  stream  are  used  first.  Thus,  it  is  straightforward  to 
calculate  the  cost  associated  with  the  exergy  units  removed  from  a 
material  stream  in  a  process  step  provided  that  this  cost  is 
calculated  in  the  previous  steps  during  which  the  exergy  units 
currently  removed  from  the  stream  are  then  supplied  to  it. 

Specific  Exergy  Costing/ Aver  age  Cost  { SPECO/AVCO )  Approach. 
This  method  was  first  presented  by  Lazzaretto  and  Tsatsaronis 
[48,52].  Following  improvements  of  the  methodology  by  Tsatsar¬ 
onis  and  co-workers  it  includes  combination  with  fuzzy  inference 
systems  for  a  more  exact  exergoeconomic  evaluation  of  plant 
components  [54]  and  systematic  improvement  of  SPECO  metho¬ 
dology  [71].  According  to  the  approach,  fuels  and  products  are 
defined  systematically  by  registering  exergy  additions  to  and 
removals  from  each  material  and  energy  stream.  The  records  of 
cost  additions  to  and  removals  from  the  same  stream  in 
conjunction  with  the  application  of  the  LIFO  principle  are  then 
used  to  calculate  local  average  costs.  This  method  consists  of  three 
main  steps:  (1)  identification  of  exergy  streams,  (2)  definition  of 
fuel  and  product  for  each  system  component,  and  (3)  cost-balance 
equations.  This  method  has  been  largely  and  successfully  used  and 
applied  to  thermal  systems  by  the  researchers  in  the  area  of 
thermoeconomics  [4,67,71-81]. 

Modified  Productive  Structure  Analysis  {MOPSA)  approach.  This 
approach  was  first  presented  by  Kim  et  al.  [45].  The  method  is 
based  on  the  exergy  costing  approach  without  flow-stream  cost 
calculations.  Lozano  and  Valero  [33]  presented  the  cost  evaluation 
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technique  of  the  thermal  system,  which  provided  a  theoretical 
basis  for  the  exergy  costing  method.  In  the  MOPSA  method,  a  cost- 
balance  equation  is  obtained  by  assigning  a  unit  exergy  cost  to  each 
disaggregated  exergy  in  the  stream  at  any  state.  The  monetary 
evaluations  of  various  exergy  costs,  as  well  as  the  production  cost 
of  electricity  for  the  thermal  system,  are  obtained  by  solving  the  set 
of  equations  for  the  unit  of  exergoeconomic  costs.  This  method  has 
been  applied  and  compared  with  SPECO  approach  for  the 
exergoeconomic  analysis  and  optimization  of  combined  heat 
and  power  plants  [45-57,60]. 

4.2.  Calculus  methods 

Calculus  methods  are  built  on  differential  equations.  Cost  flows 
in  a  system  are  developed  in  a  link  between  optimization 
procedures  that  are  based  on  the  Lagrange  multipliers  and  they 
are  used  to  determine  marginal  costs.  The  characteristic  difficulty 
in  the  application  of  calculus  methods  to  complex  systems  such  as 
combined  heat  and  power  plants  is  the  fact  that  the  Lagrange 
multipliers  vary  from  iteration  to  iteration  when  component 
thermoeconomic  isolation  is  not  achieved  [6].  Following 
approaches  are  used  in  calculus  methods. 

4.2.2.  Thermoeconomic  functional  approach  (TFA) 

This  approach  was  first  studied  as  a  Ph.D.  thesis  by  Frangopoulos 
[22-27]  while  the  first  remarkable  application  of  it  was  the  CGAM 
problem  [38].  This  methodology  is  based  on  the  Lagrangian  method 
of  mathematical  optimization.  Its  complete  implementation 
requires  the  existence  of  a  sufficiently  accurate  simulation  of  the 
system  in  order  to  determine  the  first  order  derivatives  of  the 
objective  function.  The  method  is  based  on  the  decomposition  of  the 
system  into  subcomponents,  which  may  or  may  not  correspond  to  a 
physical  component  of  the  system.  Each  component  has  a  single 
product,  and  costs  can  be  determined  by  the  solution  of  the  cost 
balances  of  the  system  components.  The  TFA  optimization  is  based 
on  the  direct  use  of  an  (usually  nonlinear)  optimization  algorithm 
since  it  requires  the  least  effort  in  the  analysis  of  complex  systems. 
The  method  introduces  the  functional  diagram  concept  and  assigns 
only  one  function  and  one  product  to  each  component  so  that 
auxiliary  equations  are  not  needed  in  the  solution  procedure.  This 
method  does  not  give  any  information  about  the  physical  and 
economic  interrelationships  among  the  system  components.  On  the 
other  hand,  scaling  of  the  variables  and  of  the  objective  function  is 
required  to  achieve  convergence  to  the  optimum  point.  Discrepan¬ 
cies  that  exist  in  searching  convergence  to  the  global  optimum  can 
be  solved  by  searching  different  initial  points  in  the  problem.  In 
order  to  find  the  optimum  solution  of  more  complex  systems, 
Frangopoulos  formulated  Intelligent  Functional  Approach  ( IF  A )  which 
was  a  further  development  of  the  TFA  [24,25]. 

4.2.2.  Engineering  functional  analysis  (EFA) 

The  basics  of  the  EFA  theory  were  developed  by  von  Spakovsky 
and  Evans  [30].  The  first  laborious  work  of  the  method  was  the 
CGAM  problem  [39].  The  basic  decomposition  algorithm  of  the  EFA 
method  is  the  variation  on  the  algorithm  proposed  by  Frangopoulos 
[22,23].  It  is  based  on  a  multi-dimensional  version  of  Modified 
Regula-Falsi  Method  [  83  ].  In  the  EFA  approach,  any  thermoeconomic 
model  exists  on  two  levels:  as  a  system  basis  model  and  as  a  set  of 
detailed  subgroup  models  [39].  These  two  models  include  informa¬ 
tion  on  the  internal  geometry  and  material  composition  of  each 
subgroup.  The  extent  to  which  these  individual  subgroup  optimums 
are  consistent  with  the  global  system  optimum  depends  on  how  well 
the  isolation  of  the  subgroups  is  established  by  the  optimization  of 
the  system  basis  model.  The  optimization  procedures  of  two 
interrelated  models  mentioned  is  possible  with  a  certain  number  of 
iterations  which  occur  in  the  link  between  the  models  and  internal 


economy  of  the  system  presented.  The  decomposition  strategy  used 
in  calculus  methods  in  thermoeconomics  is  based  on  the  Principle  of 
Thermoeconomic  Isolation  by  Evans  [6].  The  idea  is  stated  as  follows: 
A  component  of  a  thermal  system  is  thermoeconomically  isolated 
from  the  rest  of  the  system  if  its  production  P;  and  the  unit  cost  of  the 
resources  A,  are  known  quantities  and  independent  from  the  rest  of 
the  component  variables.  Decomposition  may  only  approach  the 
global  optimum  because  of  the  degree  of  thermoeconomic  isolation 
of  the  independent  variables,  the  choice  of  subgroups  and  their 
functions,  and  the  nature  of  dependent  variables  [83]. 

4.2.3.  The  structural  theory  of  thermoeconomics  (STT) 

This  was  proposed  as  a  standard  and  common  mathematical 
formulation  for  all  thermoeconomic  methodologies  employing 
thermoeconomic  models  that  can  be  expressed  by  linear  equations 
[34].  The  Theory  of  Exergy  Cost  (TEC),  the  SPECO/AVCO  approach 
and  the  Thermoeconomic  Functional  Analysis  (TFA)  can  be  dealt 
with  the  structural  theory.  The  Last-in-First-Out  (LIFO)  approach 
can  also  be  reproduced  with  the  structural  theory. 

5.  Comparison  of  the  thermoeconomic  methodologies  on  the 
CGAM  problem 

The  CGAM  problem  was  constructed  to  examine  the  optimum 
design  parameters  of  a  simple  combined  heat  and  power  plant 
which  produces  30  MW  electricity  and  14  kg/s  of  saturated  steam 
at  20  bars  [36].  The  two  optimization  criteria  were  the  maximiza¬ 
tion  of  exergetic  efficiency  and  the  minimization  of  the  total  cost 
rate  of  the  plant.  The  installation  of  the  plant  consists  of  a  gas 
turbine  followed  by  an  air  preheater  which  uses  part  of  thermal 
energy  of  the  gases  leaving  the  turbine,  and  a  heat-recovery  steam 
generator  in  which  the  required  steam  is  produced  (see  Fig.  2) 
[2,36,53].  The  base  case  design  data  of  CGAM  problem  is  given  in 
Table  3  [2,36-39,53]. 

In  the  definition  of  CGAM  problem,  three  models  were 
considered:  the  physical  model  which  describes  the  behavior  of 
the  system;  the  thermodynamic  model  which  describes  the 
equations  of  state  used  to  calculate  the  thermodynamic  properties; 
and  the  economic  model  which  describes  the  equations  for 
calculating  the  capital  costs  of  the  components.  In  order  to  simplify 
these  models,  the  following  assumptions  were  made:  (i)  the  air  and 
the  combustion  gases  behave  as  ideal  gases  with  constant  specific 
heats,  (ii)  for  combustion  calculations,  the  fuel  is  taken  to  be 
methane  CH4,  and  (iii)  all  components  except  the  combustion 
chamber  are  adiabatic. 

5.2.  Conventional  approach 

In  the  first  study  of  the  CGAM  problem  [36],  physical, 
thermodynamic  and  economic  models  were  described  in  detail 
for  each  component  of  the  plant  (see  Fig.  2).  In  this  primary  study, 
the  optimization  problem  was  expressed  as  the  minimization  of 
the  total  cost  flow  rate  of  the  plant  CT,  subject  to  the  constraints 
imposed  by  the  physical,  thermodynamic,  and  cost  models  of  the 
installation.  That  is, 

CY  =  CfrhfLHV  +  Zc  +  Zap  +  Zee  +  Zj  +  Zhrsg  (18) 

where  cf  is  the  cost  per  unit  energy  of  the  fuel  for  the  total  plant,  rh  f 
is  the  mass  flow  rate  of  the  fuel,  LHV  is  the  lower  heating  value  of 
the  fuel,  and  Zk  is  the  total  capital  cost  associated  with  the 
subcomponents  of  the  plant,  i.e.,  compressor,  air  preheater, 
combustion  chamber,  turbine  and  heat  recovery  steam  generator. 
This  first  solution  of  the  CGAM  problem  was  obtained  through 
conventional  optimization  techniques.  The  decision  variables  of 
the  CGAM  problem  were  the  compressor  pressure  ratio  P2/Pi, 
isentropic  compressor  efficiency  rjsc ,  isentropic  turbine  efficiency 
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Net  Power 


Air 


Fig.  2.  The  flow  schematic  of  the  CGAM  problem.  C:  compressor,  AP:  air  preheater,  CC:  combustion  chamber,  T:  turbine,  Ev:  evaporator,  Ec:  economizer,  ST:  steam  tank,  HRSG: 
heat  recovery  steam  generator. 


rist ,  temperature  of  the  air  entering  the  combustion  chamber  T3, 
and  temperature  of  the  combustion  products  entering  the  turbine 
T4.  All  other  thermodynamic  variables  were  calculated  as  functions 
of  the  decision  variables.  The  mathematically  optimal  solution 
results  of  this  introductory  study  [36]  are  given  in  Table  4. 

5.2.  Exergoeconomic  methodology  approach 

In  the  second  study  of  the  sequential  CGAM  problem  papers 
[37],  exergoeconomic  techniques  were  applied  and  discussed  for 
design  optimization.  The  exergoeconomic  evaluation  of  the  system 
considered  (see  Fig.  2)  consists  of  a  detailed  exergy  analysis,  an 
economic  analysis  conducted  at  the  component  level,  calculation 
of  the  cost  of  each  stream  using  an  appropriate  exergy  costing 
method  [31],  and  a  final  evaluation  at  the  component  level 


Table  3 

Base  case  design  data  of  the  CGAM  problem  [2,36-39,53]. 


Table  4 

Optimum  values  of  the  decision  variables  and  costs  in  the  CGAM 
problem  through  mathematical  optimization  method  [36]. 


conducted  with  the  aid  of  some  relevant  exergoeconomic  variables 
[2].  The  exergy  analysis  was  performed  to  identify  the  sources  of 
the  thermodynamic  inefficiencies  and  to  facilitate  development  of 
design  changes  to  improve  the  overall  system  efficiency.  Thus,  the 
exergy  destruction  in  a  component  can  be  calculated  from  the 
exergy  balance  as 

ed  =  ZC-Z  C  (19) 

‘TOT 

where  E{  is  the  flow  rate  of  total  exergy  associated  with  the  ith 
material  or  energy  stream  at  the  inlet  (for  subscript  in)  and  outlet 
(for  subscript  out )  of  the  component.  The  exergy  values  of  fuel  and 
product  for  the  kth  plant  component,  exergy  destruction  may  be 
written  as 


Ed.i<  —  Efj<  ~  Ep,k  ~  Elj<  (20) 

where  Elk  represents  the  exergy  loss  in  the  kth  component.  In 
these  papers,  the  performance  of  an  energy  system  from  the 
thermodynamic  point  of  view  was  evaluated  in  addition  to  the 
exergy  destruction,  the  exergetic  efficiency,  and  two  exergy 
destruction  ratios  [2,31,37].  The  exergetic  efficiency  is  the  ratio 
of  product  exergy  to  the  fuel  exergy: 


(21) 


In  addition  to  the  exergy  destruction  EDk  and  the  exergetic 
efficiencies  sk,  the  exergy  destruction  ratio  yDtk  is  used  in  the 
thermodynamic  evaluation  of  a  component.  This  ratio  compares 
the  exergy  destruction  in  the  kth  component  with  the  total  fuel 
exergy  supplied  EF  totai  to  the  overall  system: 


yDj< 


&F, total 


(22) 


Alternatively,  the  exergy  destruction  rate  of  the  kth  component 
can  be  compared  to  the  total  exergy  destruction  rate  EDjtota/ : 


Ep,k 

Epjotal 


(23) 


The  exergy  loss  ratio  is  defined  similarly  to  Eq.  (22),  by 
comparing  the  exergy  loss  to  the  total  fuel  exergy  supplied  to  the 
overall  system 


Yl, total 


E-L, total 
E F, total 


(24) 
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Table  5 

Comparison  of  the  global  optimum  values  of  decision  variables  and  costs  with  the  two  workable  design  condition  optimal  values  in  the  CGAM  problem  through 
exergoeconomic  analysis  and  optimization  methodology  [37]. 


Decision  variables 

Optimum  case 

Case  14  [37] 

Case  22  [37] 

p2/p, 

8.52 

10.0 

8.0 

0  sc 

0.847 

0.850 

0.850 

Ost 

0.879 

0.875 

0.870 

t3  (K) 

914.3 

900.0 

933.0 

t4(  K) 

1492.6 

1485.0 

1500.0 

Component 

Ar  (%) 

A  £(%) 

Ar  (%) 

Ae  (%) 

Ar  (%) 

Ae  (%) 

Air  compressor 

4.08 

1.30 

9.04 

1.82 

4.18 

1.31 

Air  preheater 

22.03 

-5.64 

12.73 

-4.27 

25.28 

-5.93 

Gas  turbine 

4.56 

0.87 

1.35 

0.50 

2.88 

0.73 

HRSG 

36.71 

-15.26 

38.07 

-15.58 

37.02 

-15.30 

h/rof  ($/h) 

46.70 

47.4 

49.4 

DlZ  (S/li) 

26.70 

26.8 

27.2 

Estimated  value  of  the  objective  function  ($/h) 

1302.8 

1307.9 

1307.6 

Calculated  value  of  the  objective  function  ($/h) 

1303.4 

1308.8 

1305.3 

The  detailed  economic  analysis  of  the  CGAM  problem  presented 
in  the  second  study  [37]  was  conducted  at  the  component  level  as 
part  of  the  exergoeconomic  evaluation  and  simplified  by  assuming 
some  economic  factors  such  as  capital  recovery  factor  (/3),  annual 
system  operation  time  at  the  nominal  capacity  (r)  to  specific 
numerical  values.  Since  the  exact  formulations  of  the  cost  balances 
and  additional  auxiliary  equations  depend  on  the  exergy  costing 
method  applied  (see  Eqs.  (1)  through  (7))  and  on  the  types  of 
exergy  values  (i.e.,  total,  chemical,  physical,  thermal,  and 
mechanical  exergy),  two  separate  exergy  costing  methods  were 
applied  [37].  These  methods  were  the  traditional  average-cost 
exergy  costing  method  (AVCO)  [7,8]  and  the  specific-cost  exergy 
costing  method  (SPECO)  [48,52].  AVCO  method  was  the  simplest 
method  applied  to  CGAM  problem  among  the  others.  The  SPECO 
method  was  based  on  a  detailed  accounting  of  the  cost  at  which 
every  exergy  addition  to  a  material  stream  occurs.  Thus  in  a 
process  during  which  exergy  was  removed  from  a  stream,  the  cost 
at  which  the  removed  exergy  had  been  previously  supplied  to  the 
stream  was  calculated  and  was  charged  to  the  product  of  the 
process  [48,52,54,71].  These  two  methods  were  applied  to  the 
problem  combined  with  (i)  physical  and  chemical  exergy  values 
and  (ii)  thermal,  mechanical  and  chemical  exergy  values.  In  this 
paper,  two  major  exergoeconomic  variables  were  defined  for  the 
first  time:  the  relative  cost  difference  factor  (Ar)  and  the 
exergoeconomic  factor  (f)  which  would  be  very  useful  parameters 
in  the  evaluation  of  exergoeconomic  analysis  and  optimization  of 
thermal  systems  later  [37,45,46,51,52,54,60,73,77,80,81].  The 
relative  cost  difference  rk  for  the  kth  component  is  defined  as 


that  cost  savings  in  the  entire  system  might  be  achieved  by 
improving  the  component  efficiency  (reducing  the  exergy  destruc¬ 
tion)  even  if  the  capital  investment  for  this  component  will 
increase.  On  the  other  hand,  a  high  value  of  this  factor  suggests  a 
decrease  in  the  investment  costs  of  this  component  at  the  expense 
of  its  exergetic  efficiency  [2].  Three  optimal  cases  were  compared: 
global  optimum  case  of  the  CGAM  problem  and  two  different  initial 
workable  design  case  optimal  solutions,  namely  cases  14  and  22 
[37].  Comparison  of  the  results  for  these  three  cases  is  given  in 
Table  5.  The  exergoeconomic  analysis  and  optimization  technique 
developed  and  applied  to  CGAM  problem  in  the  second  study  of  the 
sequential  papers  show  that  the  method  is  a  valuable  tool  in  the 
exergoeconomic  analysis  of  the  complex  existing  thermal  systems 
and  it  has  been  applied  to  different  thermal  systems  in  the 
numerous  papers  [1,4,31,33,37,45,46,51,52,54,60,73,77,80,81] 
since  it  was  published. 

5.3.  Thermoeconomical  functional  analysis  approach  (TFA) 

The  application  of  the  thermoeconomical  functional  approach 
(TFA)  to  the  CGAM  problem  was  first  studied  by  Frangopoulos  [38]. 
He  applied  three  different  analysis  and  optimization  procedures  in 
his  study:  (i)  direct  use  of  a  nonlinear  programming  algorithm,  (ii) 
thermoeconomic  functional  approach,  and  (iii)  modular  simula¬ 
tion  and  optimization  of  the  system.  In  the  direct  application  of  the 
nonlinear  programming  method,  the  objective  of  the  optimization 
was  defined  as  the  minimization  of  the  total  cost  rate  (i.e.  capital 
and  operation  expenses)  of  the  system 


rk 


cP,k  —  cF,k 
cF,k 


n 

(25) 

min  F  =  y^Zr  +  c  frii  fLHVu 

r=  1 


(27) 


In  an  iterative  cost  optimization  of  a  system,  if  the  cost  of  fuel  of 
a  major  component  changes  from  one  iteration  to  the  next,  the 
objective  of  the  cost  optimization  of  the  component  should  be  to 
minimize  the  relative  cost  difference  instead  of  minimizing  the 
cost  per  exergy  unit  of  the  product  with  this  component  [2].  In 
evaluating  the  performance  of  a  component,  we  want  to  know  the 
relative  significance  of  each  category.  This  is  provided  by  the 
exergoeconomic  factor  fk  defined  for  the  kth  component  as 


Zk  +  CF.k(^D.k  +  Euc) 

The  total  cost  rate  causing  the  increase  in  the  unit  cost  from  fuel 
to  product  is  given  by  the  denominator  in  Eq.  (26).  Accordingly,  the 
exergoeconomic  factor  expresses  as  a  ratio  the  contribution  of  the 
non-exergy  related  cost  to  total  cost  increase.  A  low  value  of  the 
exergoeconomic  factor  calculated  for  a  major  component  suggests 


where  Z,  represents  amortized  capital  cost  rate  of  component  r, 
including  fixed  charges  and  maintenance.  A  specific  computer 
program  was  used  for  the  numerical  solution  of  the  optimization 
problem  consisting  of  the  main  program,  double  precision 
function,  constraints  subroutine  and  the  optimization  algorithm 
developed  by  Frangopoulos  himself.  In  the  statement  of  the 
functional  optimization  problem,  the  same  objective  function 
(Eq.  (27))  was  used.  In  the  two  methods,  the  same  independent 
variables  were  taken  into  account  as  in  the  previous  applications  of 
conventional  optimization  method  [36]  and  exergoeconomic 
optimization  [37].  In  the  solution  of  the  functional  optimization 
problem,  Lagrange  multipliers  were  used.  The  simplified  form  of 
the  Lagrangian  was  defined  for  the  CGAM  problem  as  [38] 

6 

L  =  r  —  ^rYr)  +  0.2  —  ^0.2yo.2)  +  ^3^3.0  +  ^53^5.0  (28) 

r=l 
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where 

6 

rr  =  Zr  +  X>r'.rW.r,  r  =  1,2, ...  ,6 

r' 


the  case,  and  therefore  modular  approach  was  applied  to  the 
problem  only  with  demonstrative  purposes. 

5.4.  Engineering  functional  analysis  approach  (EFA) 


The  required  conditions 

VxL(x,y,  A)  =  0,  VyL(x,y,  A)  =  0,  VxL(x,y ,  A)  =  0  (30) 

were  imposed  on  the  Lagrangian  to  obtain  the  equations 


d 


s£r'  =  ° 


r=  1 


(31) 


von  Spakovsky  constructed  a  thermoeconomic  structure  for  the 
CGAM  problem  by  applying  engineering  functional  analysis 
method  (EFA)  [39].  Methodology  in  his  paper  consisted  of  two 
distinct  parts:  an  economic  model  which  represents  the  thermo¬ 
dynamic  and  economic  behavior  of  the  system  and  an  algorithm  to 
solve  the  model  for  optimum.  The  general  form  of  objective 
function  for  CGAM  problem  in  this  paper  was  defined  as  [39] 


A 


0.2  = 


dr  Q2 

dyo.2 


A 


r 


drr 

~dyf 


A 


r' 


.r  — 


A 


r 


(32) 


minimize  ftotal  —  Cc  +  CAPH  +  Ccc  +  CT  +  CHRSG  +  r  fuel  +  if  (35) 
with  respect  to 


(33)  *  =  {*,-}  =  {xc},  i  =  1,2, ...,/  and  c  =  C,APH,CC,T,HRSG  (36) 

subject  to  system  products,  component  inputs,  junctions, 

(34)  branches,  and  physical  limits  respectively  as 


In  the  solution  of  the  nonlinear  system  presented  for  the  CGAM 
problem,  numerical  methods  were  used  and  the  minimization  of 
Lagrangian  (Eq.  (28))  with  respect  to  x  by  means  of  a  nonlinear 
programming  algorithm  that  was  found  to  be  more  efficient  than 
the  analytic  derivatives  (by  avoiding  infeasible  points)  procedure 
[38].  The  modular  approach  to  the  CGAM  problem  was  given  for 
demonstration  purpose  in  the  paper  since  the  system  presented 
through  the  problem  is  not  very  complex  and  detailed  design 
characteristics  are  not  necessary.  The  comparable  optimization 
results  for  the  nominal  values  of  parameters  in  terms  of  three 
methods  presented  [38]  are  given  in  Table  6. 

Among  three  methodologies  applied  by  Frangopoulos  in  his 
paper  as  part  of  sequential  studies  of  the  CGAM  problem  in  the 
same  issue,  the  direct  use  of  nonlinear  algorithm  was  the  simplest 
one  because  it  required  the  least  effort  in  system  analysis,  but  it  did 
not  include  any  information  about  the  internal  economy  of  the 
system  (i.e.  physical  and  economic  relationships  among  the 
subcomponents  of  the  CGAM  system)  [38].  The  CGAM  problem 
was  solved  by  TFA  successfully  since  the  distribution  of  functions 
developed  for  subcomponents  of  the  system  established  inter¬ 
relations  between  components  or  between  system  and  the 
environment  leading  to  a  functional  diagram  of  the  system 
[22,23].  However  the  capabilities  of  TFA  and  IFA  which  was  a 
further  development  of  TFA  are  better  revealed  in  more 
complicated  situations  than  CGAM  problem,  when  the  solution 
of  the  complete  optimization  problem  is  required  [22-25,38].  The 
modular  approach  has  some  significant  advantages  when  detailed 
design  of  components  is  required.  This  methodology  can  be 
combined  with  direct  nonlinear  algorithms  and  TFA  or  IFA  method; 
in  particular  decomposition  of  the  optimization  problem  is 
applied.  In  Frangopoulos  paper  for  CGAM  problem,  this  was  not 

Table  6 

Comparison  of  the  nominal  values  of  the  decision  variables  and  total  cost  of  the 
plant  for  the  CGAM  problem  through  direct  use  of  the  nonlinear  programming 
method,  thermoeconomical  functional  approach  and  modular  approach  [38]. 


Decision  variables 

Method 

Direct  use  of  the 
nonlinear  algorithm  [38] 

TFA 

Modular  approach 

p2/Pi 

8.59730 

8.59770 

8.59050 

Osc 

0.84641 

0.84650 

0.84653 

rjst 

0.87886 

0.87871 

0.87878 

T3  (I<) 

912.77 

913.14 

912.93 

t4  (K) 

1491.40 

1491.97 

1491.50 

F  ($/ year) 

1.0426  x  107 

1.0426  x  107 

1.0426  x  107 

hpj  =  Yj-Yj(Xpj)=0,  j  =  1,2  ...  ,P  (37) 

hc  =  Yj  -  Yj(xc,  Yc)  =  0,  j  =  P  +  1  . . . ,  C  and  c 

=  C,APH,CC,T,HRSG  (38) 


2 


YJm  -  ^YJ+in  =  0, 

in= 1 

2 

j  =  c  + 1.. 

. ,  L  and  m  =  1 . 

,.,L 

(39) 

E  YBnout  ~  Yj  =  0, 

j  =  L  +  1. 

. .  ,B  and  n  =  1  . , 

(40) 

out= 1 


gk(x,y) > 0,  k  =  1  ...,K  (41) 

where 

Q  =  C*(xc,  Yc),  c  —  C,APH ,  CC,  T,  HRSG  (42) 

C  fai  =  r$pel(Xcc,  YCc)  (43) 

Each  of  the  general  mathematical  functions  defined  in  Eqs.  (35) 
through  (43)  were  derived  based  on  the  equations  given  in  the 
introductory  paper  of  CGAM  problem  series  [36,39].  The  internal 
economy  of  the  system  was  determined  by  EFA  using  Lagrange’s 
method  of  undetermined  multipliers  since  with  this  approach,  the 
real  unit  or  marginal  costs  associated  with  developed  thermo¬ 
economic  structure  and  model  in  the  paper  were  calculated.  Thus 
Lagrangian  was  constructed  from  the  objective  function  and 
equality  constraints  such  that  [39] 

total  =  E^  +  -  yiK)]  -  E  XAY> 

c=  1  j=  1  j=p+ 1 

-  Vj(xc,vc)]  -  £  kj[Y]m -&»<*) 

j=n+l,m=l  \  in=l  J 

-  E  ij(±y*--Yj)  (44) 

j=L+l,n=l  \out=l  J 

The  methods  used  in  the  optimization  process  in  terms  of  EFA 
method  were  categorized  in  the  paper  as:  modular,  the  Lagrangian 
and  decomposition  approaches.  For  decomposition  approaches, 
two  different  sets  of  subgroups  were  designated  namely  Case  1  and 
Case  2  [39].  In  Case  1,  the  first  set  of  subgroups  divided  the  system 
into  two  subgroups:  the  air  compressor  as  one  and  the  remaining 
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Table  7 

Optimum  values  of  the  objective  function,  component  and  fuel  costs  and  the 
decision  variable  set  for  different  methodologies  of  optimization  through  EFA  [39]. 


Design 

Modular 

Lagrangian 

Decomposition 

Decomposition 

variables 

approach 

approach 

approach  (Case  1 ) 

approach  (Case  2) 

p2/p, 

8.59822 

8.59859 

8.49996 

8.49996 

Osc 

0.84678 

0.84675 

0.84991 

0.84991 

Ost 

912.601 

912.281 

924.764 

900.09 

T3(K) 

0.87874 

0.87900 

0.88000 

0.87991 

k(K) 

1491.566 

1491.440 

1528.970 

1480.153 

^total  (SIS) 

0.36204 

0.36204 

0.36850 

0.36336 

Ca,  (US) 

Aw  ($/s) 

0.03653 

0.03658 

0.04740 

0.03467 

0.32552 

0.32546 

0.32110 

0.32869 

components  as  the  other 


*suM=j_  „,cj  and  xsub2  =  {T3,T4,t]It}  (45) 

In  Case  2,  the  second  set  of  subgroups  divided  the  system  into 
four  groups:  the  air  compressor  as  one,  the  air  preheater, 
economizer  and  evaporator  as  another,  the  combustion  chamber 
as  third  and  the  gas  turbine  as  fourth.  Its  decision  sets  were 


Xsub  1 


Xsub2  —  {^3}  5  Xsub3  —  { ^4 }  ?  K sub4 


(46) 


A  comparison  of  some  of  the  major  results  obtained  with  each 
approach  in  terms  of  EFA  method  is  given  in  Table  7.  Optimum 
results  for  modular  and  Lagrangian  were  obtained  as  almost 
identical  because  these  two  methods  coincided  at  the  optimum 
(see  Table  7)  [39].  The  decomposition  approaches  for  both  cases 
had  some  significant  differences.  This  was  expected  since  the 
decision  variables  sets  defined  for  each  subgroup  were  not  disjoint. 
In  Case  1,  more  realistic  results  were  obtained  compared  with  Case 
2  although  four  subgroups  decomposition  solution  would  be 
expected  better  than  two  subgroups  solution.  The  analysis  and 
optimization  methodologies  applied  through  this  paper  showed 
that  EFA  methodology  is  a  very  general  one  and  therefore  very 
useful  when  applied  to  the  modeling,  optimization  and  analysis  of 
energy  systems  either  time  dependent  or  not  [39]. 


5.5.  Exergetic  cost  theory  (ECT  or  TEC)  approach 


Exergetic  cost  theory  (ECT)  [33]  formulation  of  the  optimization 
problem  for  CGAM  system  was  based  on  the  use  of  Lagrange 
multipliers  [40].  The  optimization  problem,  using  the  structural 
information  provided  for  CGAM  problem  in  the  paper  was 
formulated  in  a  compact  form  as 

minimize  77  =  tceF  +  tuZ  (47) 

subject  to  constraints 

I<dP  =  F(  I),  Ps  +  (PF)F  =  P(II),  Ps  =  datum(III)  (48) 

k=fk(r),  r=fr(  t),  (w=/<u(t)  (49) 

where  77  is  the  total  cost  of  the  plant  product,  ce  is  a  vector  which 
contains  the  unit  price  of  the  resources  entering  the  plant;  Fe  and  Ps 
are  the  vectors  containing  the  exergy  of  the  entering  resources  and 
of  the  products  leaving  the  plant  respectively;  u  is  a  unit  vector;  Z  is 
the  vector  containing  the  cost  of  the  operating  and  maintenance  of 
each  component;  I<D  is  a  diagonal  matrix  containing  the  unit 
exergy  consumption  of  each  subsystem;  P  and  F  are  vectors 
containing  the  exergy  of  the  product  and  of  exergy  resources 
consumed  (fuel)  of  each  component;  and  (PF)  is  a  matrix  which 
defines  the  productive  interconnection  of  the  components  [33,40]. 


Table  8 

Optimization  results  for  decision  variables  and  costs  in  the 
CGAM  problem  by  ECT  methodology  [40]. 


The  optimization  problem  with  constraints  was  transformed  into 
one  without  constraints  as 

minimize  L(X,  A)  =  I7(X)  +  Y)  A jGj{X )  (50) 

where  77  is  the  original  objective  function  of  the  CGAM  problem,  X  is  a 
vector  containing  the  variables  of  the  problem  and  Xj  is  the  Lagrange 
multiplier  associated  with  the  constraint  Gj.  It  was  demonstrated  in 
the  paper  that  Lagrange  multipliers  associated  with  the  constraints 
had  a  profound  economic  significance  [40].  Optimization  results  of 
the  CGAM  problem  by  ECT  is  given  in  Table  8  [40]. 

The  structural  information  obtained  at  the  optimum  through 
the  ECT  method  is  very  useful  when  carrying  out  a  thermo- 
economic  analysis  of  the  plant.  Since  CGAM  problem  is  not  a 
complex  type  of  problem,  the  structural  information  developed 
through  the  analysis  and  optimization  can  be  inferred  from 
conventional  thermodynamic  analysis.  In  complex  energy  systems 
it  would  be  very  difficult  [40]. 

5.6.  Modified  productive  structure  analysis  (MOPSA)  approach 

Kim  et  al.  [45]  proposed  and  applied  a  thermoeconomic  method 
called  modified  productive  structure  analysis  (MOPSA)  to  the  CGAM 
system  and  emphasized  how  the  cost  structure  of  the  cogeneration 
system  was  affected  by  the  chosen  level  of  aggregation  which 
specifies  the  subsystems.  Kwak  et  al.  [60]  investigated  the  cost 
structure  of  this  predefined  cogeneration  system  with  the  method 
proposed  by  Kim  et  al.  [45]  and  he  rewrote  the  exergy  balance 
equations  by  reflecting  the  exergy  losses  due  to  heat  transfer.  Thus, 
assigning  a  unit  exergy  cost  to  every  decomposed  exergy  stream  in 
the  cogeneration  system,  non-adiabatic  exergetic  cost-balance 
equations  were  written  as  [60] 


Co  + 


Cp  +  Tof^S,.- 

Jnlet  outlet 


Cs 


(51) 


The  exergetic  cost-balance  equation  given  in  Eq.  (51)  yields  the 
productive  structure  of  the  thermal  system  [45,51]  that  was 
developed  by  Lozano  and  Valero  [33].  Kwak  et  al.  [60]  described 
three  levels  of  aggregation  cases  in  their  paper.  In  the  first  case  (Case 

I)  which  was  the  simplest  one,  there  was  no  internal  parameter  to 
predict  the  unit  cost  of  products,  and  it  yielded  the  lowest  unit  cost  of 
electricity  and  the  highest  unit  cost  of  steam.  The  second  case  (Case 

II)  was  the  next  simplest  one.  However,  in  this  case  there  was  only 
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Table  9 

The  unit  costs  and  the  cost  flow  rates  of  the  electricity  and  steam  of  the  CGAM 
system  by  applying  MOPSA  approach  in  three  different  levels  of  aggregation  cases 
[60]. 


Cost  structure 

Case  I 

Case  II 

Case  III 

Unit  cost  of  electricity  ($/GJ) 

8.12 

8.58 

8.46 

Cost  flow  rate  ($/h) 

(879.96) 

(926.64) 

(913.68) 

Unit  cost  of  steam  ($/GJ) 

11.62 

10.55 

10.83 

Cost  flow  rate  ($/h) 

(533.24) 

(484.14) 

(496.99) 

one  internal  parameter,  the  unit  cost  of  negentropy  (i.e.,  negative 
entropy  which  is  defined  as  the  entropy  that  the  system  exports  to 
maintain  its  own  entropy  level  low),  and  with  this  level  of 
aggregation  case  the  highest  unit  cost  of  electricity  and  the  lowest 
unit  cost  of  steam  were  obtained.  Case  III  represented  the  highest 
level  of  disaggregation.  It  had  three  internal  parameters;  the  unit 
costs  of  thermal  and  mechanical  exergy  and  negentropy.  This 
yielded  average  results  between  Case  I  and  Case  II.  In  Table  9, 
comparative  results  of  three  aggregation  levels  for  the  CGAM 
problem  are  given.  The  production  costs  for  the  CGAM  listed  in 
Table  9  satisfies  the  following  overall  cost-balance  equation  for  the 
system.  Note  that  the  unit  cost  of  electricity  and  steam  cannot  be 
evaluated  from  this  equation  [45,51,60]. 

■CHE  „  r'R  r-  X  '  7  1  a  j  r-  1  „  r  r^\ 

Ex  m)  +  ExCm  +  /  =  WnetCw  +  \EX  |CBq  (52) 


Table  10 

Comparison  of  exergy  costing  methodologies  for  the  most  disaggregated  CGAM 
system  in  terms  of  the  unit  costs  of  electricity  and  steam. 


Methodology 

SPECO 

[37] 

Structural 
theory  (STT) 
[34] 

MOPSA 
(Case  III) 
[60] 

Exergetic 
cost  theory 
(ECT)  [33] 

Unit  cost  of  electricity  ($/GJ) 

7.80 

7.42 

8.46 

7.55 

Cost  flow  rate  ($/h) 

(842.40) 

(801.36) 

(913.68) 

(815.18) 

Unit  cost  of  steam  ($/GJ) 

10.45 

10.97 

10.83 

10.64 

Cost  flow  rate  ($/h) 

(479.60) 

(503.46) 

(496.99) 

(488.06) 

In  Table  10,  various  exergy  costing  methodologies  for  the  most 
disaggregated  CGAM  system  are  compared  in  terms  of  unit  costs  of 
electricity  and  steam.  The  unit  costs  of  products  of  CGAM  system 
obtained  from  exergetic  cost  theory  (ECT)  and  SPECO  approach 
[33,37]  were  close  to  the  values  of  the  simplest  aggregation  level 
(Case  I). 

As  it  was  confirmed  in  related  works  [33,34,37,51,60],  the  unit 
cost  of  products  is  dependent  on  the  level  of  aggregation  of  the 
system  when  the  irreversibility  occurred  in  subsystems  is  given 
explicitly  in  the  cost-balance  equations.  This  situation  suggests 
that  the  cost  structure  of  the  system  is  mostly  affected  by  the 
irreversibility  represented  by  the  entropy  production  rate  at  each 
component,  and  that  cost  structure  of  the  system  is  interrelated  by 
the  irreversibility  occurred  at  each  component  [60]. 


Table  11 

Some  major  works  on  the  thermoeconomic  analysis  and  optimization  of  combined  heat  and  power  production  systems. 


Author(s) 

Thermoeconomic  methodology 

Type  of  system 

Tsatsaronis  et  al.  [1  ] 

LIFO 

Simple  gas  turbine  cogeneration 

Tsatsaronis  and  Ho-Park  [4] 

SPECO/AVCO 

Simple  gas  turbine  cogeneration 

Frangopoulos  [22] 

TFA 

Complex  thermal  system 

Frangopoulos  [23] 

TFA 

Steam  cycle  cogeneration 

Frangopoulos  [24,25] 

IFA 

Steam  cycle  cogeneration 

Frangopoulos  [26,27] 

TFA  and  optimization 

Steam  cycle  cogeneration 

von  Spakovsky  and  Curti  [29] 

EFA 

Cogeneration/heat  pump  facility 

von  Spakovsky  and  Evans  [30] 

EFA 

Combined  cycle  cogeneration 

Tsatsaronis  and  Winhold  [31] 

EEA 

Combined  cycle  cogeneration 

Lozano  and  Valero  [33] 

TEC 

Simple  thermal  power  plant 

Erlach  et  al.  [34] 

STT 

Combined  cycle  cogeneration 

Valero  et  al.  [36] 

Conventional  solution 

CGAM  problem 

Tsatsaronis  and  Pisa  [37] 

SPECO/AVCO 

CGAM  problem 

Frangopoulos  [38] 

TFA 

CGAM  problem 

von  Spakovsky  [39] 

EFA 

CGAM  problem 

Valero  et  al.  [40] 

TEC 

CGAM  problem 

Lazzaretto  et  al.  [41  ] 

Thermoeconomic  diagnosis  (TD) 

TADEUS  problem 

Verda  [42] 

TD 

TADEUS  problem 

Verda  and  Borchiellini  [43] 

STT  and  TD 

TADEUS  problem 

Hua  et  al.  [44] 

EEA 

Combined  cycle  cogeneration 

Kim  et  al.  [45] 

MOPSA 

Gas  turbine  cogeneration  system 

Cerqueira  and  Nebra  [46] 

SPECO/TFA/TEC/TECD 

Simple  gas  turbine  cogeneration 

Manolas  et  al.  [47] 

TFA  and  genetic  algorithm  (GA) 

Industrial  cogeneration  system 

Lazzaretto  and  Tsatsaronis  [48] 

EEA 

Combined  cycle  cogeneration 

Kwon  et  al.  [51  ] 

SPECO/MOPSA 

Gas  turbine  cogeneration  system 

Lazzaretto  and  Tsatsaronis  [52] 

SPECO/EEA 

Combined  cycle  cogeneration 

Cziesla  and  Tsatsaronis  [54] 

IEA— fuzzy  logic 

Thermal  power  plant 

Kwak  et  al.  [57] 

EEA 

Thermal  power  plant 

Kwak  et  al.  [60] 

MOPSA 

CGAM  problem 

Valdes  et  al.  [59] 

TFA  and  GA 

Combined  cycle  gas  turbine  plant 

Silveira  and  Tuna  [62,63] 

EEA  and  optimization 

Combined  cycle  cogeneration 

Vieira  et  al.  [65,75] 

STT 

CGAM  problem 

Zhang  et  al.  [67] 

STT 

Coal  fired  thermal  power  plant 

Erdil  [68] 

EEA  and  optimization 

Combined  cycle  cogeneration 

Mazur  [69] 

EEA— fuzzy  logic 

Combined  cycle  cogeneration 

Bonnet  et  al.  [70] 

EEA 

Micro-cogeneration 

Colpan  and  Yesin  [73] 

SPECO 

Combined  cycle  cogeneration 

Li  et  al.  [76] 

TFA  and  optimization 

Dist.  multi-generation  systems 

Ozgener  et  al.  [77] 

SPECO 

Geothermal  district  heating 

Zaleta-Aquilar  et  al.  [78] 

EEA 

Steam  turbine 

Unver  and  Kilic  [79] 

TEA 

Combined  cycle  cogeneration 

Abusoglu  and  Kanoglu  [80,81] 

SPECO 

Diesel  engine  cogeneration 
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6.  Exergoeconomic  analysis  and  optimization  of  combined 
heat  and  power  production  (CHPP) 

Combined  heat  and  power  production  (or  cogeneration)  is  the 
simultaneous  production  of  power  and  usable  heat  from  the  same 
fuel  or  energy  source  [77,79-81,84].  Cogeneration  typically 
improves  the  overall  efficiency  of  energy  use,  leading  to  reduced 
fuel  use  and  pollutant  emissions,  and  is  an  example  of  distributed 
generation  [84,85].  Depending  on  the  purpose,  exergoeconomic 
studies  on  combined  heat  and  power  production  in  literature  can  be 
grouped  in  two  categories:  (1)  if  the  purpose  is  the  allocation  and 
determination  of  the  production  costs  and/or  cost  of  losses  of  the 
plant,  the  study  is  focused  on  the  monetary  flow  rate  accounting 
through  components  in  the  plant  [4,31,48,57,66,67,70,73,77- 
81,84].  (2)  If  the  purpose  is  the  optimization  of  a  power  system, 
the  study  is  focused  on  the  selection  of  the  best  conditions  for 
operating  the  system  [23-29,36-45,47,50,51,53-55,59-63,64, 
68,75].  Thermoeconomic  concepts  and  analysis  and  optimization 
of  combined  heat  and  power  production  have  been  studied  in  a 
satisfactory  amount  of  works.  Among  these  publications,  some 
specific  and  remarkable  ones  are  given  together  with  the  thermo¬ 
economic  method  used  and  type  of  system  studied  (Table  1 1 ).  In  the 
published  papers  presented  in  Table  11,  several  thermoeconomic 
methods  were  applied  to  obtain  the  unit  costs  of  the  products  for 
combined  heat  and  power  systems.  Application  of  different  methods 
to  a  given  system  normally  yields  different  values  of  product  cost 
because  the  assumptions  made  in  the  formulations  of  cost-balance 
equations  are  different  in  each  method.  Besides,  the  cost  structure  of 
combined  heat  and  power  systems  is  affected  by  the  chosen  level  of 
aggregation  that  specifies  the  subsystems.  The  level  of  aggregation 
depends  on  the  number  of  parameters  performed  in  the  determina¬ 
tion  of  the  unit  costs  of  products  [2].  Note  that  if  the  irreversibilities 
in  system  components  are  not  considered  the  unit  cost  of  the 
products  does  not  depend  on  the  chosen  level  of  aggregation  [60]. 

A  major  fraction  of  published  works  listed  in  Table  11  were 
performed  using  algebraic  methods  [1,4,31,33,37,45,46,51,52,54, 
60,73,77,80,81].  In  these  studies,  SPECO  and  MOPSA  approaches 
were  applied  to  combined  heat  and  power  systems  in  order  to 
reduce  the  subjectivity  of  fuel  and  product  definitions  and  cost 
partitioning.  This  is  done  by  the  help  of  an  iterative  exergoeco¬ 
nomic  performance  improvement  procedure  based  on  exergoe¬ 
conomic  variables  such  as  relative  cost  difference,  exergoeconomic 
factor  and  exergetic  efficiency  as  introduced  by  Tsatsaronis 
[1,2,4,31,33,37]. 

The  remaining  studies  were  performed  using  calculus  methods 
[22-27,29-34,38,39,44,48,59,70,76,78]  mainly  developed  by  Fran- 
gopoulos  [22-27]  and  von  Spakovsky  [29,30,39].  As  mentioned 
before,  these  methods  are  based  on  Lagrange  multipliers  and  they 
vary  from  iteration  to  iteration  when  thermoeconomic  isolation  of 
component  is  not  achieved.  This  was  the  main  difficulty  of  the 
method  in  the  application.  After  the  introduction  of  functional 
diagram  concept  and  by  assigning  only  one  function  and  one 
product  to  each  component  by  Thermoeconomic  Functional 
Analysis  (TFA)  method  [22-27],  the  complex  nature  of  the  calculus 
methods  were  overcome.  The  Lagrange  multipliers  show  relatively 
smaller  variations  in  the  optimization  process  with  the  application 
of  second  law  in  the  TFA  method  [60]. 

In  order  to  reduce  the  complexity  in  the  optimization  of 
complex  thermal  systems,  several  decomposition  methods  based 
on  the  second  law  were  proposed  [44,54,62,63,69].  Combined  heat 
and  power  systems  have  some  large  scale  problems  due  to  their 
complicated  nonlinear  characteristics.  Thermodynamic  restric¬ 
tions  make  the  optimization  problem  difficult  to  be  solved.  Many 
researchers  agree  that  focusing  on  the  dominant  decision  variables 
is  more  important  than  manipulating  all  variables  simultaneously 
in  the  optimization  case  [60]. 


7.  Conclusions 

Thermoeconomics  assesses  the  cost  of  consumed  resources, 
money  and  system  irreversibilities  in  terms  of  the  overall 
production  processes.  Methodologies  presented  through  thermo¬ 
economics  may  help  to  point  out  how  resources  may  be  used 
effectively  for  sustainable  development.  Assessing  the  cost  of  the 
flow  streams  and  processes  in  a  complex  system  helps  to 
understand  the  process  of  cost  formation  from  the  input  resources 
to  the  final  products.  These  analyses  can  solve  problems  related  to 
complex  energy  systems  such  as  combined  heat  and  power 
production  (CHPP)  that  could  not  be  solved  by  using  conventional 
energy  analysis. 
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